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Abstract 

The Casimir interaction between two perfectly conducting, infinite, concentric cylinders is com- 
puted using a semiclassical approximation that takes into account families of classical periodic 
orbits that reflect off both cylinders. It is then compared with the exact result obtained by the 
mode-by-mode summation technique. We analyze the validity of the semiclassical approximation 
and show that it improves the results obtained through the proximity theorem. 
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I. INTRODUCTION 



The existence of an attractive force between two uncharged, perfectly conducting parallel 
plates was predicted by Casimir more than fifty years ago[|l|. Such force has recently been 
measured at the 15% precision level using cantilevers A similar force between a conduct- 
ing plane and a sphere has also been measured with progressively higher precision in the last 
years using torsion balances ||, atomic force microscopes and capacitance bridges [Q. 
As described in Ref . 0, , Casimir forces may be relevant in nanotechnology. The increasing 
experimental precision revives the interest in the theoretical calculation of Casimir forces 
and energies for different geometries (see Ref.[[7[] for a recent review of experimental and 
theoretical developments). Only for a few geometries the exact results are known: parallel 
plates m, self-energy for spheres || [| [ToJ] and cylinders [JTTJ, [12], [13| . For the force between 



a plane and a sphere the exact result is not known, and only an estimation valid when both 
are close enough is available. This estimation is based on the so called proximity theorem 
fl4| . Similar estimations have been recently obtained for the force between two spheres |15" 



The problems of two concentric spheres |16| and two concentric cylinders |17| have also been 



considered recently. For both cases expressions for the Casimir energy have been obtained 
using the Abel-Plana formalism. However a detailed numerical calculation and analysis of 
the results is still missing. 

In this paper we compute the Casimir energy for two perfectly conducting and concentric 
cylinders, using approximate semiclassical methods, and compare the results with those of 
an exact calculation based on a mode-by-mode summation method. We will show that 
the semiclassical approximation describes accurately the interaction energy between the 
cylinders far beyond the range of validity of the proximity theorem. 

Let us consider a system of conducting shells S a , and denote by A a suitable set of 
parameters describing the geometry of the system. For convenience the system is enclosed 
into a large box, whose boundary £ will eventually be removed to infinity. The Casimir 
energy can be formally defined as |18 



E c = £ (A)-£ (A ) , (1) 

where Eq(A) denotes the zero point energy of the electromagnetic field for the geometry 
under consideration, and E (A ) is the one corresponding to the reference vacuum. For 
semiclassical calculations it will be enough to take the interior of S without additional 



conductors as the reference vacuum However, for the exact calculation we will use 

Eq.(|I]) with Ao describing a system in which the shells are very far away. 

In terms of the modes of the electromagnetic field, the Casimir energy can be written as 

E c = ^EK-Wp) > ( 2 ) 
z p 

where w p are the eigenfrequencies of the electromagnetic field satisfying perfect conductor 
boundary conditions on S a , and w p are those corresponding to the reference vacuum. The 
subindex p denotes the set of quantum numbers associated to each eigenfrequency. 

For some symmetric systems the Casimir energy can be computed from the sum over 
modes, Eq.(§). The mode- by-mode summation technique introduced in Ref. [f| is based on 
the use of Cauchy's theorem to convert the sum over eigenfrequencies into a contour integral, 
and it turns out to be a very efficient method of calculation in this context. 

Alternatively, the Casimir energy can be computed from the knowledge of the density of 
electromagnetic modes p(E) inside S. Of course this is as difficult as computing the sum over 
modes. However, semiclassical estimates for p(E) allow for semiclassical approximations for 
the Casimir energy |20| . The Casimir energy can be written as 

r°° 1 

E c = - E Ap(E) dE , (3) 
Jo I 

where 

Ap(-E) = Pint,a + Pext ~ Pvac , (4) 
a 

describes the change in the density of electromagnetic modes inside the box £ when the 
conducting shells S a are introduced. The spectral density in each region r = int, ext or 
vac (see Fig.|TJ), p r (E) = Y. n ^( E ~ ^w r J, is given in terms of a sum on the respective 
eigenfrequencies ui rn . 

Given the Casimir energy, one can compute the forces acting on the shells S a by taking 
appropriate derivatives with respect to their positions. For the particular case we will 
consider in this paper, namely two infinite concentric cylinders of radii a and b with a < b, 
it is easy to see using symmetry considerations that the net force on each shell vanishes. 
However, the pressure is different from zero. For example, the pressure on the inner shell is 
given by 

Pc ~ M~da~ {b) 
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where E c is the Casimir energy for a length i. 



In Ref . [|T5| , p0| , the Casimir energy for certain systems has been evaluated semiclassically 
using the fact that in the limit h — > 0, Ap(E) can be computed as a sum over periodic 
orbits of a free particle bouncing off the surfaces S a . This approximation has been applied 
to (chaotic) non-symmetric configurations and the contribution of only one isolated periodic 
trajectory has been considered. This is the case for the force between a plane and a sphere, 
or for the force between non-concentric spheres. Here we will extend the formalism to sym- 
metric configurations (integrable cavities) in which a whole family of non isolated periodic 
trajectories contribute significantly in the semiclassical limit. 

The paper is organized as follows. In Section for those readers not acquainted with 
semiclassical methods, we present an overview of periodic orbit theory and semiclassical 
approximations. Section ITTTI is devoted to the semiclassical calculation of Casimir interaction 



between two coaxial cylindrical conductors. In Section [iy] we compute the exact Casimir 
energy for the same system, using the mode by mode summation technique. In Section 
we present a detailed comparison between the exact, the semiclassical, and the proximity 
theorem results. Section |VJ contains the main conclusions of our work. 

II. OVERVIEW OF SEMICLASSICAL METHODS AND PERIODIC ORBIT THE- 
ORY: GENERAL FORMALISM 

Periodic orbit theory relates oscillations in the quantum level density of a given Hamil- 
tonian to the periodic orbits in the corresponding classical system. It has its origin in the 
well known decomposition of the spectral density as 

p{E)=p{E) +p° sc (E)- (6) 

that has a rigorous meaning only in the semiclassical (h — > 0) regime for which the scales of 
variation of p(E) and p osc (E) decouple. 

To leading order in h the smooth term p(E) is the Thomas-Fermi contribution that takes 
into account the volume of accessible classical phase space at energy E. For non relativistic 
particles confined in cavities with reflecting walls, corrections to p{E) of higher order in h 
were first given by Weyl [^T|. Later Balian and Bloch |22| derived the generalized Weyl 



formula for cavities with arbitrary smooth convex boundaries and for general classes of 
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boundary conditions. We refer the reader to the book of Baltes and Hilf |23) for many 
details and examples on the Weyl's expansion. 

The deviations from the smooth behavior are given by the oscillating contribution p osc (E), 
which might be written as a sum over classical periodic orbits. The foundations of periodic 
orbit theory (POT) are closely related to the semiclassical quantization of integrable systems 
introduced by Bohr and Sommerfeld. Using a multiple reflection expansion for the time 
independent Green function and employing the principle of stationary phase, Balian and 
Bloch |24 derived in 1972 a trace formula for cavities with ideal reflecting walls giving explicit 



results for spherical cavities. Later Berry and Tabor |25j obtained exactly the same formula 
for integrable cavities using the torus quantization rule and the Poisson summation formula 
to express the spectral density p(E) as a multiple integral and evaluating the integrals by 
the saddle point method. 

The Bohr-Sommerfeld quantization rule fails for quantum Hamiltonian systems whose 
classical counterpart displays chaotic motion. In 1971 Guztwiller presented his famous trace 
formula (see Ref. |2(J for a review) obtained from Feynman's path integrals and based on 
the semiclassical expansion for the single particle propagator. The trace formula obtained 
by Gutzwiller is only applicable if all the involved periodic orbits are isolated in phase 
space, being for this reason particularly successful in its applications to chaotic systems. 
Nevertheless it fails for systems which have degenerate families of non-isolated periodic 
orbits, such as integrable systems. 

POT has been extensively used during the last twenty years to obtain semiclassical esti- 
mates for energy spectra of quantum (many-body) systems like atoms and nuclei, with vast 
applications in nuclear and atomic physics. More recently in condensed matter physics, sev- 
eral transport properties of mesoscopic systems like metal clusters, metallic grains, quantum 



dots and microstructures have been studied employing the semiclassical tools [[27j1 . 

In recent works the semiclassical Gutzwiller trace formula was exploited to evaluate the 
leading order Casimir energy for some "chaotic" configurations of pair of conductors, like 
two conducting spheres separated a small distance apart or an open shell and a sphere 



|T5| , |20f . Although the authors of these works comment on the impossibility to employ 
Gutzwiller method to evaluate the Casimir energy when non-isolated periodic orbits are 
present, they did not discuss the possible use of the Balian-Bloch (BB) or Berry- Tabor 
(BT) trace formulae. As we will show below, they can be properly adapted to evaluate the 
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Casimir effect for integrable configurations of conductors (i.e. spherical shells, cylinders, 
parallel plates, etc.) 



Our approach will be to employ a modified version of the BT trace formula [p5 |, that 
provides the appropriate path to evaluate semiclassically the Casimir effect for the case of 
two coaxial cylindrical conductors. 

We should also mention that Balian and Duplantier|2B| obtained the distribution of elec- 
tromagnetic modes inside a conducting cavity using a multiple scattering expansion for the 
electromagnetic Green functions, that leads to a decomposition of the spectral density in a 
smooth and an oscillatory term that can be written as a sum over closed polygons. In a 
subsequent work they tested the convergence of their method evaluating the Casimir energy 



for a perfectly conducting sphere [IS 



p osc will turn to be the main ingredient in the semiclassical evaluation of the Casimir 



energy Eq. (H) and can be written, to leading order in h, as a sum [26 



P° SC (E) = ME) sm(S t (E)/h + , (7) 

n t 

running over periodic orbits labeled by t, where St is the classical action of the periodic 
orbit t whose period is T = dS t /dE. For photons of energy E, the dispersion relation is 
E = p c, being p the momentum and c the velocity of light. The classical action along a 
periodic orbit t is then 

S t (E) = fp.dq = pL t = ^L t , (8) 

with L t the length of the periodic orbit. 

The phase m is the so-called Maslov index that counts the number of caustics along the 
periodic orbit. One should distinguish between Dirichlet and Neumann boundary conditions. 
For Dirichlet boundary conditions, an additional phase ir should be taken into account in 
\i t for each bounce of the trajectory on the walls of the cavity. 

The exponent v and the amplitudes A t depend on the type of periodic orbit. For in- 
tegrable systems with d degrees of freedom, the periodic orbits are non-isolated and form 
(d— l)-parameter families, with the resulting exponent v — (d +l)/2. In chaotic systems the 
orbits are isolated and their contribution is semiclassicaly smaller with v — 1, irrespective 
of the dimensionality. 



6 



III. SEMICLASSICAL RESULTS FOR TWO COAXIAL CYLINDRICAL CON- 
DUCTORS 



We consider two infinite and concentric cylindrical conducting shells G\ and C*2, of radii 
a and b (with a < b). Let us call p\2 the spectral density for photons confined in the region 
between C\ and C2 and p°f c the corresponding leading order oscillating contribution. To 
obtain p°| c one can derive a la BT a trace formula for photons inside the conducting shells, 
starting from the Hamiltonian H = E = p c, employing Poisson summation formula and 
using the torus quantization rule ||29|| . However we find it easier to proceed in two steps. 
Firstly, if one knows the oscillating contribution p osc (E m ) for non-relativistic particles of 
mass m and energy E m = p m 2 /2m confined in a given cavity, the oscillating contribution 
for photons confined in the same cavity is obtained from p osc (E m ) after replacing E m — > 
p/c and p m /m — > c. This is not restricted to the oscillating contribution, and the same 
replacement could be done in the complete spectral density p. Secondly, for cavities with 
axial symmetry, and in particular for the geometry under consideration, the periodic orbits 
(PO) are contained in planes perpendicular to the z axis (we are considering cylinders of 
infinite length). Therefore if one knows the oscillating contribution to the spectral density in 
the bidimensional annular region between two disks, Pq c , it its straightforward to obtain P§| 
that, for a given length £, 

k k 



p™%E = he k) = — / p%\^-k^)dk z . (9) 

a C7i jo ^ — ^ 



We stress that the above equation is valid for cavities with axial symmetry, but with generic 
transverse sections. 

The oscillating contribution for non-relativistic particles confined in bidimensional annu- 
lar regions, pQ C (E m ) has been previously derived (see e.g. Richter's book in Ref. [27|] and 
Appendix A). In this configuration there are two types of PO (see Fig. |2|): those which do 
not touch the inner disk (type-I), and those which do hit it (type-II). 

The type-I orbits are polygons that may be uniquely labeled by two integers (v , w) where 
v is the number of vertices (or sides) and w the winding number around the center. For 
w = 1 one has ordinary regular polygons, and for w > 1 star-like polygons, except if 
v and w are coprimes. In this situation the label (v, w) = (n j, n k) describes the n th 
repetition of the primitive orbit with j vertices and winding number k. Type-I orbits must 
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fulfill v > v(w) = Int[um I arccos(a/6)], that is the minimum value of v depends on w and 
the ratio a/b between the radii of the inner and outer cylinders, respectively. This is a 
geometrical restriction due to the fact that cos(tcw/v) > a/b. The lengths L vw of these 
orbits are 

L vw = 2v bsm(nw/v ) . (10) 

Type-II trajectories are labeled by (v,w), where now w is the number of turns around the 
inner disk in order to come back to the initial point after v bounces on the outer circle of 
radius b. We have the same restriction v > v(w) as for type-I trajectories and the length is 
given by 

L^ym — 2 V b 



^1+ (£f -2|oob(W«))- (11) 

Following the steps described above, we write p°2° = p°^j + Pyf.n as a sum of two terms 
that take into account the contributions of type-I and type-II PO, respectively. We find (see 
Appendix A) 

P^/(JE0 = EE/« A OT J5 sin L w ) . (13) 



The +(— ) in Eq.(|12D corresponds to Dirichlet (Neumann) boundary conditions (be). On 
the other hand, there is no additional phase and therefore no difference between Dirichlet 
and Neumann be in the type-II contribution, Eq.(^). The coefficient f vw is given by f vw = 2, 
except for the self-retracing type-II orbit, that has w = and v = 2, for which /20 — 1- In 



Eq.fllBD we have defined 



A — 



* ^1 — — cos(7rw/t>)^ y- cos(nw/v) — ^— ^ j . (14) 



Having obtained the explicit expression for the oscillating contribution to the spectral 
density in the region between two coaxial cylinders we can use this expression to estimate 
the semiclassical contribution to the Casimir interaction in that region. For this purpose, 
we need to rewrite Eq.(^) for the present configuration of conductors. The internal region 
is now formed by the two conducting shells. We take the outer cylindrical shell C2 as the 
shell that limits the internal and external regions. Therefore we write, 
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Pint = Pi + Pl2 , (15) 

where p\ is the density of electromagnetic modes inside the cylindrical shell C\. p\2 has been 
already defined and its oscillating contribution computed in Eqs.fll^) and ( |l3|) . 

In order to employ the semiclassical decomposition of the spectral densities in each region 
r, Eq.(|]), we derive the smooth contributions p r {E). Starting from the Weyl expansion for 



the non-relativistic case [£l], |2J and following again the steps described above, we obtain 

ME ^2^r v ^sW7w s ' +omh (16) 

where V r and S r are the volume and the surface area of each region r. The — (+) sign in the 
surface term corresponds to Dirichlet (Neumann) be. 

Therefore replacing Eqs.@, (|l|) and fllBf) in Eq.(f|) we obtain, to leading order in h, 



Ap sem (E) = Ap(E) + Ap° sc (E) = p? c (E) + p&{E) + p^(E) - p°^{E) , (17) 

where the label sem is employed to emphasize the explicit use of the semiclassical decom- 
position, Eq.(^|). The cancellation, to leading order in h, of the net smooth contribution 
Ap(E) is a consequence of the identity Vi + Vn + V ex t — V mc = 

For the geometry under study, and when considering the total electromagnetic (e.m.) 
contribution, it is satisfied that Ap(E) = to next-to-leading order in h. This is due to the 
opposite sign of the surface term in Eq . (|i6|) for the Dirichlet and Neumann be. 

As already mentioned, in order to compute the Casimir energy Eq.(|3]) one has to remove 
to infinity the boundary S that enclose the reference vacuum. In this limit p° e s ^ t {E) — > 
and p° s a c c (E) — > since no PO survive in the limit S — ► oo. Taking this into account and in 
order to compare the semiclassical results with the exact calculations, we define in analogy 
to Eq.© 

E sem _ I r E Ap sem (E) dE =- r E (pT c (E) + pT 2 c {E)) dE = E° sc + E™ c . (18) 
2 Jo 2 Jo 

It is a well known result that in the case of cylindrical geometries the Casimir energy for 
the electromagnetic case (e.m.) can be computed as the sum of the Dirichlet and Neumann 
scalar contributions. Thus, one has 
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Tpsem rposc , rposc rposc , rposc , rposc , rposc f i c\\ 

h = e 1 + h l2 = e xd + e in + e X2D + e 12jN , (iyj 

where the subscript D(N) stands for Dirichlet (Neumann) be. 

We now turn our attention to p° sc (E), the oscillating contribution for the density of 
modes inside C\. We derive p° sc (E) adapting the previous result Eq . (|l"2"D . The PO are the 
same type-I polygons labeled by two integers (v, w). Being a the radius of Cj, the length of 
these orbits is now L vw = 2v asm(irw/v), the shortest one corresponds to the self-retracing 
diametrical orbit that has v — 2, w — 1. For a cylindrical shell like C\, it is satisfied that 
v > 2w, that is v = 2 w. In other words, the cylindrical cavity is the annular one with a 
vanishing inner radius. Therefore, 

r?m - £ £ *- ^ £ e si „ (£ L„„ ± + , (20) 

w>l v>2w \ ) 

where again the +(— ) sign corresponds to he = D(bc = N). The prefactor is #21 = 1 for the 
self- retracing orbit and g vw = 2 otherwise. 
Taking into account that 

Eft = Eft + Eft =\j™E (pft(E) + pft(E)) dE , (21) 

we perform the energy integration employing Eq.(|20"D and an exponential cutoff. It is 
straightforward to show that 

poo f E ~ V 7T 7T \ f°° / E ~ V TV 7T \ 

Jo ^ [%c L ^ ± - + 2) dE = fii/o 6XP ( " A E) ^ [he L ^ ± ^ + 2) dE 
= ±2 (^-) sin(7r v/2) , (22) 



and in conclusion 



E osc = E osc +E osc= Q (33) 



Thus in the semiclassical approximation, the e.m. Casimir energy for an isolated cylinder 
vanishes. This is also the case to lowest order in the multiple scattering expansion |19| . 

To compute E° 2 C = E° 2 C D + E° 2 C N we should take into account, in principle, both type-I 
and type-II PO for each scalar contribution, D or N. However, as before, the contributions 
from type-I PO cancel out and only type-II PO contribute to the semiclassical e.m. Casimir 
interaction. We compute E° 2 C noting that Pi 2 c n ls the same f° r Dirichlet and Neumann be 
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and therefore 

E n= 2 \j™ E P°i2,n{E)dE. (24) 

Replacing Pil c n from Eq.(|l3|) and performing again the energy integration using an ex- 
ponential cutoff, we obtain the main result of this section. Namely, 

E sem = E osc = a W J J /» 1 ^ («» ^ W ) > ( 25 ) 

4 77 a w>0 v>v V 

where we have defined a = b/a, and 

</(a- cos(t™/w)) (a cos(vrw/w) - 1) 

N(a,v,«j) =-5= 2 : —72 • (26) 

(1 + cr — 2a cos{7Tw/v)) 

It is interesting to rewrite Eq.fl2"5|) by separating the terms with w = from those with 

w > 1. The orbits with w = and t> > 1 are the diametrical orbit and its repetitions, whose 

lengths are L v0 = 2 v (b — a). Therefore we write 

E sem = E sen + ^ = _ Q l/2 £ £ J_ ^ ^ w) (2?) 

Z " fl io>l u>-0 U 

In order to obtain E^™ , we calculate from Eq . (^6|) 

1 

(a - 1 



N(a,v,0) = — ^, (28) 



and using the fact that J2 v >i v 4 = C( 4 ) = ?r 4 / 90 we finally obtain, 



tt =° 360 a 2 (a - l) 3 " 360 (6 - a) 3 ' 

This expression, which is valid for arbitrary values of a > 1, diverges as a — > 1. One 

might wonder whether the sum of the contributions with w > 1 cancels and/or contributes, 

in any way, to this divergence. To investigate this problem it is interesting to note that, in 

such a limit, it is possible to derive an analytic expression for the sum of these contributions. 

Due to the fact that for a ~ 1, 



v = Int[wn/ arccos(a/6)] ~ ° != » w , (30) 

P (a - 1) 



we can approximate cosfaw/v) ~ 1 — (n w/v) /2. Thus, from Eq.(p26 

M« - hv,w) = . (31) 

(a — I) 6 („,2 I Om) 2 a x ~ 

V V U ~+~ (a- 1)2 
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Since v » 1 we can convert the sum J2v>i> ^ N(a ~ l,v,w) in Eg. ([27]) into an integral. 
Performing the integration and summing over w > 1 we obtain 

. hc C(3) 



- M ' " ^ " + £>(« -1), (32) 



4 7r 3 a 2 a 



which implies that remains finite as a — > 1. Thus, in this limit the semiclassical result 
is divergent and completely dominated by the w = contribution. As we will see later, this 
dominance extends over a rather large range of values of a. 

Having shown that for a ~ 1 one has E sem ~ E!^™ it is interesting to compare Eq.(p9|) 



with what results from the application of the proximity theorem ||14|| . The Casimir energy 
between two parallel plates of area A separated a small distance 5 is 

^ h c it 2 A , . 

For the two coaxial cylinders, when (b - a) C a ~ 6, the proximity approximation for the 
Casimir energy, E P , is obtained taking 5 = (b — a) and A — > 2 ti£ a ~ 2 n tb in Eq. (|33|) . 
The result is 



h C 7T 3 



for A = 2 7i£ a 



E P ={ 360 ( fe " a ) 3 ' (34) 

The proximity theorem does not suggest any particular choice for the area A. This is irrel- 
evant in the limit a ~ 1. 

The semiclassical result Eq . (p9|) can be interpreted as a proximity approximation with a 
given choice for the area: the geometric mean between the areas of both cylinders. Although 
it has been shown that to compute a proximity force for a "chaotic" geometry, the geometric 



mean radius should be taken ||30| , such demonstration does not apply to the present case. 
However, the semiclassical approximation provides a justification for it. As we will see 
in Section V, this choice reproduces the exact pressure with an error less than 10% for 
1 < a < 4, improving the proximity result beyond its expected range of validity. 

IV. EXACT CASIMIR ENERGY 

In this section we will compute the exact Casimir energy for the coaxial cylinders using 
the mode by mode summation method of Ref. 0. 
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The zero point energy in Eq.fl2|) is ill defined, because both series are divergent. In order 
to define it properly we introduce a cutoff function as follows 

EeM = I J2( e ~ aWpw p ~ e-°**w p ) ■ (35) 
1 v 

The exact Casimir energy E ex is the limit of E ex (a) as a — > 0. For simplicity we choose an 
exponential cutoff, although the explicit form is not relevant. In our definition we take w. p 
as the eigenfrequencies of the electromagnetic field satisfying perfect conductor boundary 
conditions at r = a, r = b and r = R, and w p are those corresponding to the boundary 
conditions at r — Ri, r — R2 and r = R, in the limit R > R2 > Ri ^> a > b. R\, R2 and 
R are the parameters that define the reference vacuum (see Eq. (ED)). R corresponds in this 
case to the radius of the surface S in Fig.|l|. 

In cylindrical coordinates, the eigenfunctions are of the form 

h nkz (t, r, d, z) = e { - iw ^ t+ik > z+ine) R n (\r) , (36) 

where the function R n is a combination of Bessel functions satisfying the perfect conductor 
boundary conditions. These boundary conditions define the possible values of the constant 
A. The eigenfrequencies are w n k z = C\J k 2 z + A 2 . 

The frequencies of the TE modes are defined by the following relations: 

J n (Aa) = r < a 

J n (Xa)N n (Xb) - J n (\b)N n (\a) = a<r <b 

J n (Xb)N n (XR) - J n (XR)N n (Xb) = b<r <R. (37) 

For later use we introduce the notation F^ E for the product of these three relations 

F^ E (z,a,b,R) = J n {za)[J n {za)N n {zb) - J n {zb)N n {za)][J n {zb)N n {zR) - J n {zR)N n {zb)} . 

(38) 

The frequencies of the TM modes involve derivatives of the Bessel functions: 
J' n (X nm a) =0 r <a 

J' n (X nm a)N' n (X nm b) - J' n (X nm b)N' n {X nm a) =0 a<r <b 
J' n {X nm b)N' n {X nm R) - J' n {XnmR)K{X nm b) = b<r<R. (39) 
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As before, we introduce the notation 



^ M (z,a,b,R) = J' n {za)[J' n {za)N' n {zb) - J' n {zb)N' n {za)][J' n {zb)N' n {zR) - J' n {zR)N' n {zb)} . 

(40) 



The set of quantum numbers p in Eq. ( |55| ) is given by (n, m, k z ), where m denotes the 
different solutions \ nm of both F^ E (z,a,b, R) = and F^ M (z,a,b, R) = . As we are 
considering infinite cylinders, we can replace the sum over k z by an integral. The result for 
the Casimir energy for a length i is 

E ex (a) = ^ jT ^ £ (^Le--v^L _ ^ + ^ e -«V^) , (41) 

where we have denoted with A nm the solutions of the equations F^ E (z, Ri,R 2 ,R) = and 
F™(z,R 1 ,R 2 ,R) = 0. 

Using Cauchy's theorem it follows that 

hlc izze ~°T^ !{z) ^ e ~'"> (42) 

where f(z) is an analytic function within the closed contour C, with simple zeros at x\, x 2 , 
within C. Following Refs. ||, |l(J, we use this result to replace the sum over m in Eq.(f4l"l) 
by a contour integral 

E ^ - § £ f ? / ^v^-^s >- « ■ <«) 

where 

p em (zab)= lim m 

n 2cyl[ ' ' J " ft^oo F^(z, ifc, R 2 , R)F^ M (z, R u R 2 , R) ' l44j 
It proves to be convenient to compute the difference between the energy of the system of 
two concentric cylinders and the energy of two isolated cylinders of radii a and b 

E 12 (a) = E c {a) - E 1 (a, a) - E x {a, b) . (45) 

The energy of isolated cylinders has been computed previously using other regularizations 



TT] , [1211 . The formal expression with the cutoff can be obtained adapting the two cylinders 
case above. It is given by 
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where 



Fn lcyli Z i a ) 



Hm FV TE {z,a,R)FV™{z,a,R) 
r^oo F W TE ^ Ru R)Fi l)TM (z, R u R) 



(47) 



and 



F^ TE (z,a,R) = J n (za)[J n (za)N n {zR) - J n {zR)N n (za)} 
F^™(z,a,R) = J' n {za)[J' n {za)N' n {zR) - J' n {zR)N' n {za)) 



(48) 



Therefore, the explicit form for £"12 is 
£hc r°° dk 



^ = S £ f ? L tz^e-'^i * F. u (* a, 6) 



where 



n 12 



>,a,&) 



F*¥ cyl (z,a,b) 



F E ¥,M,a)F E ^,,(z,b) 



(49) 



(50) 



To proceed we must choose a contour for the integration in the complex plane. In order 



to compute E ex (a), Ei(a, a) and Ei(a,b) separately, an adequate contour is [10fl a circular 
segment Cr and two straight line segments forming an angle <fi and it — <p with respect to the 
imaginary axis (see Fig. |3|). The nonzero angle <fr is needed to show that the contribution of 
Cr vanishes in the limit r — > 00 when o > 0. For the rest of the contour, it can be shown 
that the divergences in E ex (a) are cancelled out by those of Ei(a, a) and Ei(a, b). Therefore 
in order to compute £12(0") we can set <fi = and o = 0. The contour integral reduces to an 
integral on the imaginary axis. We find 

the 



E 



12 



2na' 



12 



(iy, l,a) 



(51) 



where we have used that F n i 2 (iy/a, a, b) = F n u(iy, l,a). 

Using the asymptotic expansions of Bessel functions we obtain, on the imaginary axis 



31 



F™ yl (iy,l,a) 
F n 12 (iy, l,a) 



lQa 2 y 4 I n (y)I' n (y)[I n (y)K n (ay) - I n (ay)K n 
[I' n (y)K' n (ay) - I' n {ay)K' n {y)]K n {ocy)K' n {ay) 
-4y 2 I n (y)I f n (y)K n (y)K' n (y) 



I n (y)K n (ay) 
I n (ay)K n (y) 



I' n {ay)K' n {y) 



(52) 
(53) 
(54) 
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Note that F n n{iy-, is a real function. Therefore, the integral over y in Eq. ( [ST]) is 
restricted to y > k z . We rewrite Eq.fl5"T|) as 

£hc f°° f°° i d 

E12 — ~7T~ 2 ~ 2 I dkzY. dyJk 2 z -y 2 — ln F n 12(^,1, a) • (55) 
ztt a Jo „ J k z dy 



After some trivial steps we obtain 
the 



the f°° 
E12 = 4^2" H J dy y\nF n 12 (iy,l,a) 

ffir ( roo roo \ 

N dyy\nF 12 (iy, l,a) + 222J o dy y\nF n 12 (iy,l,a) J . (56) 



Ana? 

The exact energy for the two concentric cylinders is the sum of E\% an d the Casimir energies 
for single cylinders of radii a and b Jl2|, |I3| . The final result is 



E ex = E 12 - 0.01356 (- + -)£hc. (57) 

In the limit a ~ 1, E\ 2 dominates the Casimir energy and can be evaluated analytically. 
Using the uniform expansion for the Bessel functions we obtain, to leading order in (a — 1) 

F n 12 (iy, 1, a) « [l - e- 2 ^"- 1 )] 2 , (58) 

where 

z 2 



Hz) = 1 + — -j=== • (59) 
1 + V 1 + 2 

It is convenient to expand 



00 — 2nh(z)(a— l)k 

ln[l - e -^Hz)(*-i)] = - . (60) 

k=i k 

Inserting this expression into Eq. (^) one can compute first the sum over n, then the integral 
over y and finally the sum over k. The contribution of the n = term in Eq.(^) is negligible 
for a ~ 1. The final result is 



As expected, in this limit we obtain the proximity approximation |14[ for the Casimir energy 
(see Eq.(||)). 



16 



V. NUMERICAL RESULTS 



In this section we present the numerical results corresponding to the semiclassical and 
exact calculations described in the previous sections. 

Before comparing the results, it is worth to note that , as described in Section III, the 
semiclassical e.m. energy for an isolated cylinder vanishes (Eq. (p3[)), while the exact energy 
for a cylinder of radius a is —0.01356 I he/ 'a 2 (Eq.(|57|)). Therefore, in principle one could 
compare the semiclassical Casimir energy with the full exact energy given in Eq. fl5"T|) , or with 
the interaction energy Eyi- This leads to some ambiguities, although the second possibility 
seems more appropriate because only trajectories bouncing off between both cylinders give 
a non vanishing contribution to the semiclassical energy. 

In Fig. |] we display the exact results (dashed line) together with those corresponding 
to the semiclassical calculation (full line) as a function of the ratio of the radii of the two 
cylinders, a = b/a. For convenience we have defined the dimensionless energies as 

„2 rp „2 rpsem 

a &i2 sem _ a tj 

612 -~Yc~' e -~% c ~r- (62) 

As all energies are defined up to a constant, it is more meaningful to compare the exact 
and semiclassical pressures on a given cylinder. Specifically, we compare the pressure on the 
inner cylinder due to the presence of the outer cylinder, 

Pl2 = ~ w^T ' (63) 

with 

i a rpsem 

We also display in Fig. | the exact and semiclassical results for the dimensionless pressures 

27rfl4 sem 27Fa4 sem 

P12 = -j— Pn , P = -j— V ■ (65) 

As expected, the discrepancies between the results increases with a. One should remark 
that, as already mentioned, the comparison between the exact and semiclassical energies is 
subjected to some ambiguities that are not present in the case of the pressures. Therefore, 
assuming a typical experimental precision of the order of 10 %, such discrepancies would 
be observable only for values a > 4 if the pressure would be measured (see Fig.| lower 
panel). Thus, it is fair to conclude that, within a typical experimental error, the semiclassical 
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expressions lead to very good results for values of a in the rather wide range 1 < a < 4. It 
should be stressed that, within that range, the Casimir interaction energies and pressures 
decrease by several orders of magnitude. 
The full pressure on the inner cylinder is 

1 dE prr fie 
p ex = -TT-o^ = Pia - 0.01356— . 66 

For a ~ 1 the attraction of the outer cylinder dominates and the full pressure is positive. 
The inner cylinder tends to expand. However, when a is large, the self-pressure (second 
term in the above equation) is bigger and the full pressure becomes negative. The crossover 
takes place at a ~ 3.15. It is remarkable that the semiclassical approximation is accurate 
beyond this critical value. 

Given the good results obtained within the semiclassical approximation it is interesting to 
consider them in greater detail. We have seen that in the case of the w = contribution an 
explicit analytical expression can be obtained for arbitrary values of a. One might wonder 
how well the exact results can be described by this expression. Numerical calculation shows 
that if one would plot such contribution in Fig. [| they would be completely indistinguishable 
from those corresponding to the full semiclassical calculation. In fact, even for a = 10 the 
contributions with w > 1 represents less than 3% of the total value of the energy or pressure. 
Thus, it is clear that the semiclassical results are completely dominated by the w = 
contribution which, in turn, is dominated by the primitive self-retracing orbit (v = 1) and 
its first repetitions (v = 2, 3). We have also mentioned that the expression corresponding to 
the w = contribution, Eq. (|29|) , reduces exactly to that given by the proximity theorem 
in the limit a — > 1. However, as explained at the end of Sec. Ill, it is not clear how to 
extrapolate the use of such theorem away from that limit. The basic problem is that there 
is an ambiguity on which area should be considered for the differential surfaces. In fact, 
one might either take the area given by the radius of the inner cylinder, the one given by 
the radius of the external cylinder or any combination in between (e.g. average radius, 
geometric mean radius, etc). The situation is illustrated in Fig. [| where our semiclassical 
results (full line) are compared with the extrapolation of the proximity theorem using the 
two extreme alternatives: area given by the inner cylinder (dashed line) and that given by 
the outer cylinder (dotted line). Obviously, any other sensible choice would lie between 
these two curves. Although for values of a very close to one this ambiguity is completely 
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harmless, already at such small values a ~ 1.12 it implies an uncertainty of the order of 10 
%. The present semiclassical calculation removes completely such ambiguity indicating that 
indeed the geometric mean radius is the right choice. 



VI. CONCLUSIONS 

In this article we have computed the Casimir interaction energy for two infinite con- 
centric cylinders exactly and using an approximate semiclassical method. To perform the 
semiclassical calculation, we have extended the method introduced in Ref.|20|, valid for 



non-symmetric configurations with isolated periodic trajectories, to the case of symmetric 
configurations in which families of non isolated periodic orbits contribute significantly in the 
semiclassical limit. Technically, this involves the use of the Balian-Bloch or Berry- Tabor 
trace formulae instead of the Gutzwiller formula. As for the exact calculation, we have used 
the mode- by-mode summation method of Ref. 0, but using a cutoff instead of zeta function 
regularization. 

The final result for the semiclassical Casimir energy (Eq.fl2"o])), is dominated by the self- 
retracing periodic orbit (w = 0) and its repetitions. Moreover, it can be interpreted as a 
proximity approximation with an effective area given by the geometric mean of the areas of 
both cylinders. While this choice for the area was previously derived for proximity forces in 
non-symmetric configurations |15|, the semiclassical calculation provides a justification 
for the case of an integrable cavity. We have found that, surprisingly, the semiclassical 
result describes accurately the Casimir pressure beyond the range of validity of the proximity 
theorem. Indeed, while this theorem is expected to be valid when b — a <C a, the semiclassical 
result reproduces the exact pressure within a 10 % up to b = 4a. 

As a side point, we remark that the semiclassical approximation fails to reproduce the 
electromagnetic Casimir energy for an isolated infinite cylinder. Unlike the case of a rect- 
angular parallelepiped PD"| , the Dirichlet and Neumann contributions have opposite signs. 



Taking into account that for cavities with axial symmetry, the electromagnetic Casimir en- 
ergy is the sum of Dirichlet and Neumann scalar contributions, we obtained a vanishing 
result for the semiclassical energy (other approximations also give a null result for the cylin- 



der [|19], |32fl). This shows that the semiclassical method may not work for cavities with one 



of the dimensions much larger than the others, as was suggested in Ref. [ 20|| . 



19 



The semiclassical approach described here could be applied to other integrable configura- 
tions of conductors, like two concentric spheres. Moreover, it can be extended to incorporate 
effects that could be relevant experimentally, like small surface deformations or roughness. 
We are currently investigating these problems. 
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APPENDIX A 



In this Appendix we describe with some detail the steps followed to obtain the oscillating 
contributions p°{£i and plfu, Eqs. fll2D and (|13|). 



As we have mentioned in Sec.|ITT] the oscillating contribution for non-relativistic particles 
confined in a bidimensional annular region 0, pQ C (E m ), has been previously derived (see 



e.g. Richter's book in Ref.[27|). It can been written as a sum of two terms, each one 



associated respectively to the contributions from type-I and type-II PO's. In Sec.|T| we have 
characterized these PO, whose lengths are L vw = 2v bsin(irw/v) for type-I PO, and L vw = 
'2 r h ( / j 1 + — 2 | cos(7ra/t>) J for type-II PO (see the description that precedes Eqs. 



( [TOP and (O)). Taking this into account we can write p°Q C (E ri 
with E m = h k 2 /2 m [27J and 



P&% E m) + Pqji(E„ 



= g g f- ^ ^ cos (* L m ± ^ + I) , (Al) 
12 b 2 m ( 7r\ 

^■• {Em) = S g u 2 n ww^ A - sin (* L ™ + V ■ (A2) 

been defined in Eq.([TJ|). Therefore in order to obtain the oscillating contributions 
to the density of modes for photons confined in the region we have to replace in Eqs. ( |A1| ) 
and (|A2|) , h k 2 /2 m — > p/c and (h k) /m —> c. This is trivial and leads to Po c j(k) and pQ C u(k) 
with k = E/hc. 

To derive the oscillating contribution to the spectral density in the region between the 
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two coaxial cylinders we have to perform the integral Eq.(^), that we write here in a slightly 
different form, 



PlUE = he k) = J- £ -jJ== p%Uft^W *. , (A3) 

to emphasize with the subscript po = I or II the contributions from the two types of 
PO's. PQ^ P o{\fk 2 — kl) is obtained after replacing k — > ^jk 2 — k 2 z in p°Qj{k) and in Po°n(k). 
The integrals are straightforward to perform to leading order in ft, and as a result we obtain 
Eqs.flg) and flT|). 



21 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948); V.M. Mostepanenko and N.N. 
Trunov, The Casimir effect and its applications, Clarendon, London (1997); M. Bordag, The 
Casimir effect 50 years later, World Scientific, Singapore (1999); P. Milonni, The quantum 
vacuum, Academic Press, San Diego (1994); G. Plunien, B. Muller and W. Greiner, Phys. 
Rep. 134, 87 (1986). 

[2] G. Bressi, G. Carugno, R. Onofrio and G. Ruoso, Phys. Rev. Lett. 88, 041804 (2002). 
[3] S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997). 

[4] U. Mohideen and A. Roy, Phys. Rev. Lett. 81, 4549 (1998); B.W. Harris, F. Chen and U. 

Mohideen, Phys. Rev. A 62, 052109 (2000). 
[5] H.B. Chan et al, Science 291, 1941 (2001). 
[6] H.B. Chan et al, Phys. Rev. Lett. 87, 1801 (2001). 

[7] M. Bordag, U. Mohideen and V.M. Mostepanenko, Phys. Rep. 353, 1 (2001). 
[8] T.H. Boyer, Phys. Rev. 174, 1764 (1968). 

[9] V.V. Nesterenko and I.G. Pirozhenko, Phys. Rev. D57, 1284 (1997). 
[10] M.E. Bowers and C.R. Hagen, Phys. Rev. D 59, 02007 (1999). 
[11] L.L. DeRaad Jr. and K. Milton, Ann. Phys. (N.Y.)136, 229 (1981). 
[12] K.A. Milton, A.V. Nesterenko and V.V. Nesterenko, Phys. Rev. D 59, 105009 (1999). 
[13] P. Gosdzinsky and A. Romeo, Phys. Lett. B 441, 265 (1998). 

[14] B.V. Derjaguin and LI. Abriksova, Sov. Phys. JETP 3, 819 (1957); B.V. Derjaguin, Sci. Am. 
203, 47 (1960). 

[15] M. Schaden and L. Spruch, Phys. Rev. Lett. 84, 459 (2000). 
[16] A. A. Saharian, Phys. Rev. D63, 125007 (2001) 
[17] A. A. Saharian, ICTP preprint, IC/2000/14. 
[18] See Plunien et al, in Ref[l]. 

[19] R. Balian and B. Duplantier, Ann. Phys. 112, 165 (1978) 

[20] M. Schaden and L. Spruch, Phys. Rev. A 58, 935 (1998). 

[21] H. Weyl, Nach. Akad. Wiss. GSttingen 110, (1911). 

[22] R. Balian and C. Bloch, Ann. Phys. 60, 401 (1970); 63, 592 (1971). 

[23] H. P. Baltes and E. R. Hilf: Spectra of finite systems ( B - I Wissenschaftsverlag, Mannheim, 



22 



1976). 

[24] R. Balian and C. Bloch, Ann. Phys. 69, 76 (1972). 

[25] M. V. Berry and M. Tabor, Proc. R. Soc. London, Ser. A. 349, 101 (1976). 
[26] M. C. Gutzwiller in Chaos in Classical and Quantum Mechanics (Springer- Verlag, New York, 
1990). 

[27] For a review see for example: M. Brack and R. K. Bhaduri in Semiclassical Physics (Addison- 
Wesley Publishing Company, Masachussets, 1997); K. Richter in Semiclassical Theory of Meso- 
scopic Quantum Systems (Springer, Berlin, 2000). 

[28] R. Balian and B. Duplantier, Ann. Phys. 104, 300 (1977). 

[29] M. J. Sanchez, unpublished. 

[30] J. Blocki, J. Randrup, W.J. Swiatecki and F. Tsang, Ann. Phys.105, 427 (1977). 

[31] For simplicity, to compute the limit R,Ri,B,2 — * oo we assumed R2/R1 = ot. The final result 

for the Casimir energy does not depend on this particular choice. 
[32] G.J. Maclay, H. Fearn, P. W. Milonni fcuant-ph/0105002 



23 




FIG. 1: Conducting shells S a , not necessarily connected, limiting internal and external regions, 
with the space cut-off £ limiting the external region. 
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FIG. 2: Two periodic orbits in a bidimensional annular region. The type-I orbit (v = 4, w = 1) does 
not touch the inner disk (long-dashed line). The type-II orbit (v = 8,w = 1) hit it (short-dashed 
line). 
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FIG. 4: Dimensionless Casimir interaction energy (upper panel) and pressure (lower panel), defined 
in Eqs.(^) and (|65[), as a function of a = b/a. In both panels the dashed line corresponds to the 
exact result and the full line to the semiclassical result. 
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FIG. 5: Dimensionless Casimir interaction energy in the semiclassical approximation (full line) 
as compared to the result obtained using two naive ways of determining the relevant area which 
enters in the proximity theorem: area of the inner cylinder (dashed line) ; area of the outer cylinder 
(dotted line). 
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